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Abstract 



The thermal self-energy of an electron in a static uniform magnetic field B is cal- 
', culated to first order in the fine structure constant a and to all orders in eB. We 

use two methods, one based on the Furry picture and another based on Schwinger's 
proper-time method. As external states we consider relativistic Landau levels with 
special emphasis on the lowest Landau level. In the high-temperature limit we derive 
self-consistent dispersion relations for particle and hole excitations, showing the chiral 
asymmetry caused by the external field. For weak fields, earlier results on the ground- 
state energy and the anomalous magnetic moment are discussed and compared with 
the present analysis. In the strong-field limit the appearance of a field-independent 
imaginary part of the self-energy, related to Landau damping in the e^e~ plasma, is 
pointed out. 
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1 Introduction 



It has become increasingly important to understand the quantum field theory of elementary 
processes in the presence of both a thermal heat bath and strong background fields, in 
particular in connection with applications to astrophysical and cosmological models. A 
natural strategy for making progress in this area is to generalize existing background field 
calculations to finite temperature and vice versa. The first obvious quantity to study is 
the effective potential, which governs the thermodynamics of the system, and there are 
several papers on this topic [1-7]. Next, we are interested in more detailed kinematical 
issues, which are determined by the propagation of the elementary excitations in a plasma. 
Dispersion relations for electrons have been studied in great detail in strong background 
fields but without any heat bath [8-14], and at high temperature but in the absence of 
external fields [15-19]. 

In this paper we perform a detailed study of the fermion self-energy in QED, taking 
into account both the effects from a magnetic background field and a thermal heat bath. In 
particular we consider the limits of strong fields in the lowest Landau level and weak fields at 
high temperature. There are two basic methods for doing perturbative quantum field theory 
calculations in a constant magnetic background field. One is the Furry picture, where 
the fermion propagator is constructed by an explicit sum over the solutions to the Dirac 
equation [^. The other one is the Schwinger proper-time method, where the propagator is 
expressed directly in terms of operators without any reference to an explicit representation 
of the states [^ . The different methods are suitable for different calculations and we have 
used them both. Typically, the Furry picture is convenient in the strong field limit where 
only a few Landau levels contribute, while Schwinger's method is particularly useful for 
weak fields. 

We have organized the paper as follows. The two basic methods, the Furry picture 
and Schwinger's proper-time method, are introduced in Sections ^ and ^. The weak field 
limit is studied in Section ^ using both methods. It is particularly interesting to find the 
dispersion relation at high temperature for weak fields, as we do in Section ^ since it 
is then consistent to not only calculate first-order corrections but to solve the dispersion 
relation self-consistently. As an application of the weak field expansion we study the 
anomalous magnetic moment in Section |^ and compare it with other calculations. The 
strong field limit, where only the lowest Landau level contributes, is investigated in Section 
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||. Discussions of the results follow in Section |, and the appendices contain some technical 
details. 



2 Self-energy corrections of the lowest Landau level 
in the Furry picture 

We consider Dirac fermions with charge — e in the presence of an external field described 
by the vector potential = (0, 0, Bx, 0), corresponding to a static uniform magnetic field 
in the negative 2;-direction. Using static energy solutions to the Dirac equation + ej^ — 
m)^^^^ = 0, we may represent the second quantized fermion field in the Furry picture 
as 



vl/ = ^[6.vl.(+)(x,t) + 4vI/W(x,t) 



(2.1) 



where k denotes a complete set of quantum numbers and (6, c?) are the standard annihilation 
operators for particles and antiparticles p2| . The fermion propagator, including the effects 
of some distribution of particles, can then be constructed explicitly as the expectation 
value (see 0] for more details) 



iS{x,x) = (T [*(x)*(x)l) 



(2.2) 



Evaluating the expectation values of the creation and annihilation operators we may sepa- 
rate into vacuum and thermal (actually more generally due to some arbitrary distribution 
of fermions) contributions, S{x', x) = S^a_c{x', x) + S^'^{x', x). We find 

Y: \e{t' - t)¥^\x'wi'-\x) - Q{t - t')^i-\x')¥-\x) 



iS^'f'ix^x) 



, (2.3) 
(2.4) 



where E^^ is the energy eigenvalue of "^l^^- For fermions, the antiparticle distribution is 
determined by fp{ko) = 1 — fp{—ko). We shall here consider only the case of thermal 
equilibrium: 
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(2.5) 



where f3 is the inverse temperature and /i is the chemical potential determined by the 
charge density of electrons and positrons of the system. For the explicit form of the wave- 
functions and of the propagator we have summarized the relevant expressions in Appendix 



To lowest order in the electromagnetic coupling there are two possible contributions to 
the electron self-energy, the one-particle irreducible self-energy diagram and the tadpole. 
In Appendix |C| we show that the tadpole gives no contribution to the self-energy in a 
neutral background. Therefore, in configuration space the self-energy is to one-loop order 
given by 

- zS(a;', x) = {-lef-f^iD^'^ix' - x)iS{x\ x)-f^ . (2.6) 
In a general covariant gauge, parametrized by C,, the photon propagator is given by [23-25] 



—I 



-27rS{q')fB{qo) 



+ le 

where the photon distribution function in the case of thermal equilibrium is 



(2.7) 



exp[/5|go|] - 1] 

and the operator d/dq^ is not supposed to act on the go in fsilo)- We have written 
[d'^q] = d"-q/{2TT)"' as a short-hand notation for the integration measure. The expectation 
value, AE, of the self-energy in a state described by the wave function '^^^n,py,pS-''^^ 
defined by the expression 

.J, /ci^xdV^g,^,^(xOS(x-,x)vI^g,^,^(x) 

where / dydzdt is a normalization factor from the continuous spectrum in p^, p^, since 
/ d'^x'^^^{x)^ ^{x) = J dydzdt, cf. Eq. (|A.16| ). In the gauge we use, the translational 
invariance in the plane perpendicular to the B field is reflected by the energy degeneracy 
for different values of Py. Explicit calculations show indeed that AE(^.n,py(yPz) is independent 
of Py in the lowest Landau level, so we consider for simplicity only py = here. We formally 
show in Appendix |^ that the self-energy is independent of the gauge-fixing parameter ^ 
on the tree-level mass shell. This is explicitly verified in the high-field limit in Section |8.2| . 
We shall therefore use the Feynman gauge ^ = for the photon propagator. From the 
inverse propagator to one-loop order, the effective Dirac equation is obtained as 

(i^-m)^(x,t) =y"dVS(a;,x')^(a;') , (2.10) 

where D^^ = — ieA^. Let us consider the value of the energy E, appearing in the wave 
functions for particles in Eq. ( A. 3 ) as exp[iEt], as a free parameter not to be constrained 
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by the tree-level Dirac equation. Then the expectation value of Eq. ( |2.1CI| ) leads to the 
dispersion relation 

E = EM + ^E{p,) , (2.11) 



where E^iPz) = + + '^^Bn, and AE is to be calculated on the tree-level mass shell. 
It is tempting to try to solve Eq. ( |2.11| ) self-consistently in E, by substituting En ^ E 



in the phase factor in \l/ before taking the expectation value, and calculate AE with 
such an arbitrary E. However, this is inconsistent for several reasons. First, not only 
E but the actual form of the wave functions should then also be solved self-consistently. 
Secondly, it turns out that the self-energy is in general only gauge-invariant on the tree- 
level mass shell and a self-consistent solution would become gauge-dependent. In |^ a 
fermion propagator equivalent to the one used here, but written explicitly in terms of 7 
matrices, was derived. Using this propagator and the general ansatz for the wave functions 
in Eqs. ( |A.3| ) and (|A.10|) , it is possible to calculate the general matrix structure of the self- 
energy operator, that acts on the space-independent spinor m^, in the Furry picture. Then 
one could proceed as in the following sections, where instead a generalization of Schwinger's 
proper-time method has been employed, and solve the self-consistent dispersion relation in 
a limit, such as the high-temperature limit, where the off-shell gauge dependence may be 
neglected. 

Here we shall confine ourselves to the calculation of the energy shift for an electron 
in the lowest relativistic Landau level. Suppressing the C = l) "^ = 0, Py = subscripts 
and the B and pz dependence, we separate the self-energy into its contributions from the 
vacuum, thermal fermions, and thermal photons 

AE = AE,^ + AEX^+AE^ , (2.12) 



respectively. In [|T^ the vacuum part was calculated for p^ = 0, using the methods adop- 
ted here. Introducing an arbitrary momentum parallel to the external field does not alter 
anything in the vacuum sector, since the relativistic invariance is unbroken in the z and 
t directions. The self-energy is then a function of only Eq ~ p^ = and eB. We have 
confirmed this by explicit calculations, of which we only quote the result here. 



A_£/vac = d,s du exp(— m us ) 

Eci 2n Jo Jo 

(2.13) 



2eB 



1 + se' 



-2eBus 



1 + 2eBu{l - s) - e-2e-B«s ^ 



The second term on the right-hand side comes from the ordinary mass renormahzation to 
make the vacuum part of the self-energy vanish for i? = 0. 

Using the thermal electron propagator given in Appendix 0, we have 



AE' 



-I— - ^ , , I exp[z£'o(t - t ) - 2^42; - 2;)J/o;o(a;)/o;o(a; ) 



Eq J dy' dz' dt 
X / [dS] 



,4 ^ exp[— ■ (x — x')] 



(f + ie 



00 



xJ2 M^o] t'^^?/] l^^^] exp[-iko{t - t) + iky{y - y) + ik^{z - z)] 

.2 2 ,2 



X {m^/„;fc,,(a;)/„;^(a;') + - Eoko+p^k^ In-i-ky{x)In-i;ky{x')} 



[2.U) 



The t, y and z integrals are trivially performed, and produce (5-functions that are used to 
perform the corresponding q integrals. We now express In-ky in terms of the explicit form 
in Eq. ( |A.5|) , and use the fact that Hq{x) = 1. Then the variables x, x',k^ = and ky are 
shifted and rescaled. Fourier-transforming the (5-function 

1 1"°^ 

5{kl-m'^ -kl-2eBn) = — dsexp[-is{kl - - kl - 2eBn)] , (2.15) 

2tc J-oo 



and using the Feynman prescription to write 

a + ie Jo 

with a = {ko- E^f - {k^ - p^f - 2eB{kl + k^) + ie, we find 

?p,tj- _ _^ / dk^fpiko) I ds du dxdx' dk^dkydk^ 

I J —CO Jo J 



(2.16) 



A r.P a • 2eB 



^-00 



X expS^-is[kl - - kl] + iu[{ko - E^f - (k, - p^)^ - 2eB{kl + kyf + ie]^ 



X exp 



m 



-^i^' + ^'')-kyi^ + ^')+^kAx-x')-2kl 

00 

+ (m' - Eoko + p,k,)X\ ^ —Hn{x)Hn{x') 



n=0 



(2.17) 



where we have shifted the summation in the In~i,n-i term, and defined A = e*^*^^*. Let us 
now use the integral representation of Hermite polynomials |^ : 

2 



exp[-a;72]iy„(a;) 



dv exp[— f /2 + ivx]{—ivy 



2.18 
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and similarly for Hn{x') in terms of an integral over v'. We may then identify an exponential 
function 

f:^t^^Ip^ = eM-vv'X] . (2.19) 

We now have some Gaussian integrals that may be performed in the order v,x,v', and 
finally x'. Let us introduce a factor of convergence exp[— 2e-Be {k^ + ky)], where e — > 0+. 
We may now also change the order and perform the k^, ky integrals, while keeping the s, u 
parameter integrals. The final result then reads 

AE^f = -i — — / dkofpiko) / ds du dk^exp{-is[kl - - kl]} 
vr^ Eq J-oo J-oo Jo J 

X exp{tu[{ko - Eo) -{k,-p,) + te]} ^ , „ ^ . . . • 

(2.20) 

The Gaussian integral over kz could easily be performed, but we prefer to keep it in order 
to make the s,u integrals simpler when considering the weak field limit. 

The thermal photon contribution is obtained in a way analogous to the case of thermal 
electrons. The final result reads 



^ roc poo poo p 

- / dk^fsiko) I ds I du I c/fc^ exp{— isf/cQ — /c^]} 

J—oo J ~oo Jo J 



xexp{^MPo + ^o) -m - (k, + p,) + ^^]} ^ ^ 2eE(e - zg) - e'^^eB" ' 

(2.21) 



The two equations ( |2.20| ) and ( p.21| ) serve as starting points for the weak field expansion 
in Section ^. 

3 Schwinger's proper-time method 



In 1951 Schwinger used the analogy with a quantum mechanical problem to find an 
explicit expression for the electron propagator in an external electromagnetic field without 
first finding the solutions to the Dirac equation in the field. Later, mostly in the 70 's, this 
method was used to calculate a number of different properties of electrons in external fields. 
The real-time finite-temperature propagator can be constructed in a simple way from the 
zero-temperature propagator. Therefore, we shall in this section use Schwinger's method 
to calculate the electron self-energy in a magnetic field at finite temperature, in order to 



6 



compare with the corresponding calculations using the Furry picture, and to generalize to 
the full self-energy operator. We use standard real-time rules and as long as we are only 
interested in the real part of the self-energy to one-loop it is sufficient to calculate the 
(ll)-part of the propagator. The fact that standard Thermo Field Dynamics works also in 
a magnetic field should be obvious, and it was discussed in 0, |^. (For the imaginary part 
we can use the rules in pH].) The propagator then is 



iS{x', x") = {x'\ 



fl— m + ie 



fl— m + ie f[— m — ie 



\x } 



(3.22) 



where /f(Po) is defined in Eq. ( [A.18| ) and f[ = 7^(p^ + ey4^),as usual. We shall not re- 
peat Schwinger's calculation, but just give the result for the zero temperature part in our 
notation 



[x 



fl— m + ie 



\x") 



(An) 



{x',x") 



°o ds 



e ^^"^^ exp 



— IS 



eaF 



+ m — ie 



X exp[— -(x' — x")eF coth.eFs{x' — x")] 
X |^(eFcotheFs + eF){x' - x") + m| 
0(x',x") [[d'p]e-'P^'''-^"\SUp) , 



where we have suppressed Lorentz indices (F = F J^) and used the notation 



(3.23) 



(x\ x") 



exp 



exp[— L(s) 



dx>' ( + -F^,{x - x"Y 



2 1,2 

a — 



sin(eas) sinh(e6s) 



ah 



BE 



(3.24) 



The phase factor 0(x', x") is the only part of Eq. ( p.23| ) that depends exphcitly on the gauge, 
and that part factors out in a natural way when we compute the self-energy. Restricting 
the field to be purely magnetic in the negative z-direction we find in momentum space 



ds — 

coseijs 



exp 



ts 



P\\ 



2 tanc-Bs 2 



eBs 



— m + te 



1P\\ 



-ieBsa, 



COS eBs 



+ m 



(3.25) 



where for general four- vectors a and b, a ■ foy = ao^o ~ CLzbz and a-b± = axbx + dyby. In order 
to calculate the thermal part we choose a gauge where Aq = and OqA^ = 0. This is very 
natural for a static magnetic field (for a further discussion of the gauge dependence of the 
thermal propagator in a background field, see 0). The (f){x', x") can then also be factorized 
out in front of the thermal propagator. At finite temperature we see from Eq. ( p.22 ) that 
Eq. (|3.25| ) should be replaced by 

^^vac(p) - /f(Po) (^«5vac(p) - «5vac(p)) • (3-26) 

The real combination that occurs in the thermal part of Eq. (|3.26D is finally obtained by 



extending the s-integral in Eq. (|3.25|) from — oo to oo. In the integrand of Eq. ( p.25| ) there 



are poles and essential singularities on the real s axis. They have to be avoided by taking 
the integration contour in the lower half-plane for positive s (see e.g. for a discussion 
of this contour); to get a real quantity for the thermal part it has, therefore, to go in 
the lower half-plane also for negative s. These poles are similar to the ones responsible 
for the non-perturbative pair production in an external electric field or the de Haas- 
van Alphen oscillations in a degenerate electron gas Q]. When we consider the weak field 
limit in Section ^ we do not encounter these poles to leading order. 
The self-energy may be represented in momentum space as 

- iS(x', x") = -i(j){x', x") J [dV]e"'^^"'"""^S(p) , (3.27) 

where the only breaking of translational invariance is in the phase factor (f){x', x"). There are 
two contributions to the real part of the thermal self-energy at one loop. One from thermal 
photons and one from thermal electrons. Let us start with the photon contribution. The 
electron contribution is completely analogous. The thermal part of the photon propagator 
may be represented in the Feynman gauge by 

D^\k) = -g^^fsiko) - = -g'^fsiko) H dtexp[zte - \t\ e] , (3.28) 
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and the corresponding contribution to the self-energy is 



/roo 
[d^^fBiko) / dtexp[ite - \t\e] 
J —oo 



X I ds exp 
/o 

AeBsa, ( 



is [{p — k)t 



tan eBs 



eBs 



{p — k)\ — + ie 



-ieBscr, 



X ■ 



COS eBs 



,l{p-k)\\ - 



COS eBs 



lip - k)± + m\j^ 



le 



(2vr) 



X exp 



\dkQ\j B{kQ)dt ds 



7r 



1/2 



eBrr 



J{s + t)J i{eBt + tan eBs) 



st 



ttaneBs 



s + t eBt + tan eBs 



■p^ — m s] — se ~ \t\e 



AeBsa, 



-leBsa, 



eBt 



t e 

^ coseBsy^^^^ s + t^^"^^ cos eBs eBt + tan eBs 



7P_L + m| 7^ , 
(3.29) 



where we have followed very closely, the only essential difference being that the ko- 
integral cannot be carried out explicitly. The T,^ above is a not an operator but a com- 
plicated function of the parameter p^. In order to obtain the expectation value of the 
energy shift in a given state one would have to multiply it with explicit wave functions and 
integrate over p^. There is, however, a clever way to rewrite it as an operator (thus 
replacing p^ with gauge-invariant operators 11^) which, when acting on the tree level wave 
functions, has simple properties. Noticing that 0(x',x") only depends on x'j_ and x", we 
write 

{x'\±\x") = [ [rfp||]e-'P("'-"")ii0(x',x") / [c/px]e'^("'-""^^S(p||,px) . (3.30) 



The key relations to be used are then 



ix' x") 



exp 



—I 



V 



eB ^ 



[x I exp 

Jp{x'-x")^ 



.tant; 
eB 

x") cos V 



p± 



[dp±]e''- 



[x I exp 



exp 



tanf 



-P± 



1P± 



V 



-i—iJii-eBa,) 
eB 



%\x") cos'^v . 



(3.31) 



After performing the 7-matrix algebra we obtain the final expression for the self-energy 
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operator 



57r^ J-oo 



[dko] 



TT 



OO nOO 

dt ds { 



1/2 



eB 



eBt + tanei^s 



fsih 



X exp 



i[tkQ + s{po - /co) 



^^2 ^ tt2 2 

—11 , — m s 



COSf 

COS eBs 
2 cosf 



-2e 



eBt 



cos e5s e5t + tan eBs 



7o(Po - ^o) 



t 



s + t 



Pzlz ) + 4m cos ei?s 



where 



tanu 



eBt taxi eBs 



(3.32) 



(3.33) 



e-Bt + tan eBs 

We notice that 11^ and do not commute but is actually multiplied by a function of 
only — eBcjz with which it does commute, so the ordering is not a problem. A similar 
analysis for thermal electrons gives 

1/2 



J-oo 



[dko] 



dt 



ds 



TT 



eB 



^iis + t)^ 



eBt + tan eBs 



X exp 

COSf 



« HPo - «;o) + skQ —p^ -n^ -m s 

^ s + 1 eB 



X- 



+ 



cos eBs 
2 cos V 



-2e 



-ieBsa, 



t 



s + t 



Pzlz ) + 4m cos eBs 



eBt 



(3.34) 



cos eBs eBt + tan eBs 

The matrix elements of Sf+g- and Sf^ depend of course on the basis in which they are 
calculated. A suitable basis that diagonalizes the eigenvalues k, = {E,n,Py,Pz) should, 
with the chiral representation of 7^, have the form (see e.g. Appendix 0): 

"^^{x) = exp[-i{Et - p^z ~ Pyy)]Vn,py{x)u{E,n,py,pz) , 

Vn,py{x) = diag(^In,py{x),In~l,py{x),In,py{x),In-l,py{x)^ , (3.35) 

where In,py{x) is defined in Appendix ^ and u{E,n,py,pz) is a Dirac spinor, independent 
of x^. With this choice of wave function, the space-time integrals of (see Eq. ( |2.9| )) 

/ d^x rf^x'^^(x)S(x, x')^^-{x') 



(3.36) 



/ dydzdt 

can be carried out directly and give a 5-function in /t — n'. There then remains a 4 x 4 
matrix that can be diagonalized with a suitable choice of the spinors u{E,n,Py,Pz)- 
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4 Weak-field expansion 



In many physical applications the magnetic field is strong enough to be important, but 
still weak enough for an expansion in eB/m? to be useful. It should, however, be noted 
that there are some fundamental difficulties with the limit of weak fields. At any finite 
field strength the eigenf unctions are Landau levels and there is not a single eigenfunctions 
that continuously goes over to a plane wave when 5^0. The Fourier coefficients of the 
function /„ p,_^(x) in Eq. (|A.5|) are given by 



dxe'''''In,py{x) = ^[^] exp 



2 / TT \i/4 



{k - ipyf p. 



2 



2eB 2eB 



2 



Hn{kJ—\ , (4.37) 



and they cannot be expanded in powers of -B. It is, therefore, not clear that there exists a 
continuous limit as i? ^ in general. In fact, we know that in a background of degenerate 
electrons there are de Haas-van Alphen oscillations that do not have a series expansion in 
5 . In standard first-order perturbation theory it is enough to calculate the expectation 
value of the perturbation in the unperturbed states. It would therefore be tempting to use 
plane waves as external states for weak fields. But, as explained above, the exact states are 
not approximately equal to plane waves and we have to use the Landau levels as a basis 
even for weak fields. 



4.1 The Furry picture 

In order to obtain the weak-field limit we expand the denominator in Eq. ( p.20|) : 

0{eBy . (4.38) 



-i2eB 



eBs' 



l + 2eB(e + iu) 



oi2eBs 



U + te 



U + It) 



In the limit i? ^ we may close the contour and use Cauchy's theorem to integrate over s. 
The simple pole from the first term in Eq. ( [4. 381 ) gives a contribution for k^ — m? — k"^ > 0. 



We may then perform also the t and k^ integrals, with the final result 

dkoe{kl-m^)fF{ko) 



m 



Pz 



(4.39) 
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This result is finite, well-behaved as pz 0, and it agrees with the thermal self-energy 
calculated in a conventional plane- wave basis. To linear order in eB we use Eq. ( 4.38|) 
together with e'^*^'^'^ = l + i2eBs + 0{eBY in Eq. ( p.20|) . In order not to get a contribution 
when closing the s contour by a semi-circle at infinity, it is necessary to rewrite terms such 
as e.g. s/{s — u + ie) = 1 + u/{s — u + ie) and treat the 1 separately. The term linear in 
eB is thus obtained from 

,2 



ieB 



{Eoko -pzkz) + 2m 



u 



u 



It 



+ {2m' - Eoko+pzkz)- 



u 



u + ie) 



(4.40) 



Integrating over s the first term in Eq. ( [4.40| ) gives a contribution proportional to 5(A;q — 
— kl). For the last two terms we proceed as before, close the s-contour and then use 
Cauchy's theorem. Performing the standard integrals over t and kz we arrive at the term 
linear in eB 



a eB 
47r Eq 







Pz 



X 



-2 



kl — m? 



Pz 



2Eoko - m^) In 



Eoko - m' + pzJkl - m? 



^Eoko -w? - PzJkl - m 



{kl - m'){2kl - SEoko + m2) + pl{kl + p^ - Eoko 



(4.41) 



{ko - Eo)' 

This expression has an ostensible singularity at ko = Eo- When expanding around this 
value we find that the possible divergences cancel, and that the integral is finite as it stands 
for Pz > 0. The limit Pz ^ is considered in Section ^ 

The thermal photon contribution is obtained in a similar manner. We may in this case 
also perform the ko integrals. In the field-free case we obtain the well-known result [29-32] 

m ttT^ 



AE^{B = 0) 



a 



Eo 3m2 



(4.42) 



Similarly we find the term linear in eB 

m2 avr eB 



AE^{B) - AE^{0) 



Eo 6 pI 



Eq 

Pz 



In 



' Eo+Pz 



E, 



O-Pz, 



(4.43) 



However, these results only apply to an external electron in the lowest Landau level its 
energy being given according to the tree-level Dirac equation. In the next subsection we 
shall consider the general case. 
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4.2 Schwinger's method 



In the weak-field limit the Landau levels get closer and closer, and in the Furry picture an 
explicit method of resumming them is always necessary. This resummation was performed 
in Section ^, and the result expanded to linear order in the external field in the preceding 
section. However, this problem does not occur when using the Schwinger proper-time 
method, since no explicit reference to the Landau levels is made. In fact, with Schwinger's 
method we have an expression in terms of operators which we can apply to any state 
we wish. We shall now use this method to study the weak-field limit of the self-energy. 
Starting from Eq. ([33^ ) we find (H^ = Ul+ pi) 

i2 ^oo /•no pna ^ 



le 



47r3/2 



[dkQ 



X exp 



dt / dsfsiko) 

-oo Jo 

>o St 



(2(s + t))3/2 



Oj 



t 



rn?s 



loiPo - ko) - —zlzPz ) - 2m - — — e^^^^'^^jZl. 



- (4.44) 

after expanding most terms to linear order in B (the terms with poles as a function of s). We 
keep the full B dependence wherever it is added linearly to 11 since 

and n need not be small. In all other places the dependence is 0{B^). The t and s integrals 
in Eq. (|4.44|) can be performed using 



dt 



dt 1 



oo (zt)3/2 t 

dt 



exp 



exp 



for 6 > 0, and 



t exp 



i{ta + -) 



i{ta + 
i{ta + 



1 /2 

2('^^ sin(2v^)e(a) 



2i- 



avr 



1 1/2 /sin(2Va&) 



2Va6 




-2i 



TT 



1/2 



COS 



(2v^) e(a) 



e(a), (4.45) 



I I dsd 
h 
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i / dse''" 
Jo 



dse'^'smibs)- 
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I I dse^''"- sm{bs)i 
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sa ( smjbs) 

\ bs 
/ sin(6s) 
1^ bs ' 



cos(6s) I - 



cos(6s) 



-IS 



2 



a + b 



a'^-b'^ ' 
-1 

1 



4 In 

2b 



a — b 



2b 



In 



a 



a + b 
b 



a + b 



a?-b^ 



(4.46) 
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The final result for the photon contribution is {k = |A;o|, ^1= 7n_|_ + ^zPz 



[dko]fBiko) 





a ~ b 








d + b 


k 











[lo{po - ko 
a — b 



2m - 1/L] 



a + b 



(4.47) 



where fll 



2IIk . 



n — 



2poko 



(4.48) 



^2 



The meaning of and 11 depends on which 7-matrix they multiply, and they should be 
replaced according to 



when multiplying Il,7^,7j,, 



m + eB(Tz , when multiplying 70, 7^ 



- eBa, 



when multiplying 1,70,72 , 
when multiplying ''ix^ly ■ 



(4.49) 
(4.50) 



These are just replacement rules that simplify the writing and not any mathematical iden- 
tities. The Equation ^.47| equals exactly the self-energy in the absence of the magnetic 
field, with a replacement 11^, except that and 11^ should be shifted by ±eBaz 

according to the rules in Eq. ( [4.49| ). We write it in this short-hand way in order to more 
easily see that it reduces to the well-known result in the limit of vanishing B. At the same 
time we notice that it is not enough to simply replace by 11^, but there are some extra 
±eBaz terms that enter. The electron contribution is calculated in a similar way: 

o2 /.oo 

[dko]fF{ko)e{kl-m') 



-13 



2^ 



X 





a — b 






i2n 


a + b 



[70^0 - 2m] + — 

n 



l + ^ln 

. 2b 



a — b 



a + b 



(4.51) 



where now k = Jkq — m? and 



2 + - 2pQkQ 



Po 
2nk 



(4.52) 



Again, it is easy to check that the zero- field limit from Eq. ( 4.51 ) agrees with standard 
results. 

One advantage with Eqs. ( [4.47| , ^3TD is that they are expressed directly in terms of 
gauge-invariant operators and they can be used to calculate the expectation value between 



14 



any Landau levels. It is particularly useful in Section |^ where we solve a self-consistent 
dispersion relation without specifying the exact form of the spinors. On the other hand, 
the results agree with the ones from the Furry picture where we have checked them. An 
example is the anomalous magnetic moment that we discuss in Section ^ 

5 High-temperature limit 

Since a few years back there has been a great interest in the high-temperature limit of gauge 
theories, mainly stimulated by the successful resummation of a consistent and infinite 



set of diagrams, which is encoded in the so-called Hard Thermal Loop (HTL) ||T^, |l| 
effective action. In the high-temperature limit it is meaningful to not only compute a 
perturbative correction to the energy spectrum, but also to solve the dispersion relation self- 
consistently. The reason is that in QED the dominating 0{eT) correction comes only from 
the one-loop diagram that we have calculated, all higher-order diagrams being suppressed 
by extra factors of (9(e^T). The leading terms are also gauge-fixing independent and 
gauge-invariant, which makes the whole procedure consistent. Therefore, it is particularly 
interesting to take the high-temperature limit (T ^ m, pz, VeB, /x) and study the effects 
of a weak magnetic field. The effective Dirac equation is 

s{po, n'hoPo - r{po, n2)7,p, - r(po, - e5a,)a - ml ^ = , (5.53) 



where as before 11^ = 11^ + p^ and 

s(Po,n2) = 

r(po,n2) = 



1 



2po \n\ 



In 



Po + 


\U\ 


Po - 


\U\ 



1 



Po 

2ini 



In 



Po + 


\U\ 


Po - 


\U\ 



(5.54) 



The temperature dependence enters only through the thermal mass Ai^ = It is 

almost possible to guess the expression in Eq. ( |5.53| ) from the standard expression for the 
HTL Dirac equation. The usual momentum p^ should be replaced with the gauge-invariant 
momentum H^, but there is an ambiguity in replacing p"^ by H^ or by ]/I ]ZI . The correct 
way follows from the lengthy calculations in Section ^ It would also be interesting to 
compare Eq. (|5.53|) with the equation of motion obtained directly from the HTL effective 
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action. This is in general difficult due to the non-local and non-linear character of the HTL 
effective action, but we have checked that it agrees up to linear order in B. 

From now on we shall take the high-temperature limit and neglect the vacuum mass, 
which simplifies the Dirac equation considerably. The matrix structure of the function 
r(po, — eBaz) can be made explicit by rewriting it as 

r(po, - eBa-z) = ^(H + o^.)r(po, - e5) + ^(U - a,)r(po, + eB) . (5.55) 



In the chiral representation Eq. ( |A.1| ) we obtain the following Dirac equation 
/ 

^ = , (5.56) 



-poS-PzT -r-^+ \ 

-r+^- -Pos + Pzr 

r_e+ 

-PoS -PzT / 



where ^± = U^Ti^y and r± = r{pQ, Il^±eB). Using the ansatz in Eq. ( p.35| ), the relations 
in Eq. ( |A.8| ), and dividing \I' into left and right 2-component spinors {u = {R, L)'^ in the 
conventions of P^), we find the decomposed Dirac equations 



-ESr 



PzTn -iV2eBnr^_i 
iV2eBnr^_i -Es„„i + p^r„_i 



-ESn + PzTn 



i\/2eBnr„_i 



-iV2eBnr^_i -Esn-i - Pz^n-i 



L 



R 







(5.57) 



where s„ = s{E,pl + eB{2n + 1)), and similarly for r„. The dispersion relations follow 
immediately by taking the determinants of Eq. ( |5.57| ) 



L : {Esn + PzTn) {Esn-i - PzTn-i) " 2eBnr'i_i = 



R: {Esn-Pzrn){Es 

n-l + Pzl^ 



n—lj 



2eBnr'^ i = . 



(5.5^ 



These relations are valid for all n > 1, but in the lowest Landau level there is only one 
non-zero component for each of L and R, so the dispersion relations reduce to 



L: Es{E,pl + eB) +p,r{E,pl + eB) =0 
R : Es{E,pl + eB) -pzr{E,pl + eB) = 



(5.59) 



Since the high-temperature correction is a consistently resummed large correction to the 
vacuum dispersion relations, we want to solve them self-consistently. In the absence of a 
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Figure 1: Dispersion relation and spectral weight for the right-handed branch in the lowest 
Landau level (n = 0). All dimensionful parameters are given in units of the thermal mass Ai. 

magnetic field it is well known that the high-temperature dispersion relation for fermions 
has two branches |^ |TB|, ^ and that spin up/down (or left /right handedness) is de- 
generate. The hole branch corresponds to a positive energy solution to the factor in the 
dispersion relation that usually only gives a negative energy solution. We have a simi- 
lar phenomenon here, and the hole solution is associated with the factor Esn + Pz^n for 
positive Pz, even though the dispersion relation does not factorize completely. There are 
almost the separate symmetries {pz ^ —Pz), {E ^ —E) and {L R), but they fail be- 



cause of the difference in the index n for the two factors in Eq. (|5.58|) , which also breaks 
the spin degeneracy. We notice on the other hand that there are symmetries under the 
combinations {L ^ R, p^ ^ ~Pz)i {L ^ R,E ^ —E) and {E ^ —E,pz ^ —Pz)- The 
number of states for given values of Pz and n is eight, corresponding to {L/R)x (par- 
ticle/hole) X (positive/negative energy), apart from the usual degeneracy in py. In the 
lowest Landau level, only right-handed particles and left-handed holes can propagate for 
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0.0 I ' ' ^ ' ' ' 1 

-2.0 -1.0 0.0 1.0 2.0 

Figure 2: Dispersion relations and spectral weights for the next-to-lowest Landau level (n = 1). 
The dotted part of uJh shows a continuation of the dispersion relation beyond the point where Uh 
picks up an imaginary part. All dimensionful parameters are given in units of the thermal mass 
M. 

positive pz, and vice versa for negative pz (see Eq. ( p.59[ )), as can be understood from the 
following argument. The energy shift for a magnetic moment /x is given by — /x ■ B. The 
electron-spin contribution to ^t points opposite to the spin, since the charge is negative. 
We have chosen B to point in the negative 2;-direction so that the energy is lowest when 
the spin points in the positive 2;-direction. For pz > the helicity should thus be positive 
in the lowest Landau level and we expect the positive chirality state R to propagate (see 
Fig. |l|). Left-handed particles have the wrong helicity to be in the lowest Landau level for 
Pz > 0. However, holes have opposite chirality-helicity relation, so a left-handed hole can 
propagate for pz > 0. Similarly, right-handed holes propagate for < as shown in Fig. |l|. 
A similar asymmetry exists in higher Landau levels in the sense that, for a right-handed 
particle with momentum pz and 11^ = + eB{2n + 1), there exists another right-handed 
particle state with momentum —pz but with 11^ = pl + eB{2n — 1) and a different energy. 
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It is clear that this asymmetry can be important for decay processes where only one chi- 
rality is produced, for instance in /3-decay. To get an asymmetry of this form, we need a 
chiral charged particle with spin (or magnetic moment). It clearly needs to have a splitting 
between spin up and down, but it is also important that it is chiral in order to separate 
right- and left-handed particles at the same time. At zero temperature there are no chiral 
charged spin 1/2 particles in the Standard Model, but at high temperature the dominant 
mass effect is the thermal one, which is chirally invariant. In this way the temperature 
effect makes the electron essentially chiral without reducing its magnetic moment. 

It has been observed that the polarization of electron and positron spins in a 

strong magnetic field generates an anisotropy of the neutrino emission in weak processes, 
and that this anisotropy could be the explanation of the high space velocities of pulsars. In 
the presence of a thermal heat bath, we have seen above that a further left-right asymmetry 
is generated for the electrons and positrons. It remains to be seen to what extent this 
asymmetry increases the asymmetry in the neutrino emission during the hot phase of the 
supernova explosion. 

In order to find out the correct final-state factors for a decay into the different branches, 
not only the spectrum is needed but also the spectral weight for each branch. That is, we 



need the wave function renormalization factor Z(k). It can be computed from [17 



auj 



(Tr [{Viuj,p,,n)^or']y' , (5.60) 



where I){uJ,Pz,n) is the matrix in Eq. ( |5.56|) and Ei corresponds to the solutions for L, R, 



particle and hole. The spectral functions should satisfy the sum rule for n > 1 

Zip + Zih + + ZRh + multiparticle states = 2 . (5.61) 



We use a definition of Z that differs from by a factor of 2. In the lowest Landau level 
there are only half as many states, so the spectral weights add up to 1. The result of a 
numerical calculation of the spectral weight as well as of the spectrum and for the lowest 
Landau level is presented in Fig. [I|. It shows the right-handed branch, but the picture is 
the same for the left-handed branch with ^ ^Pz- A similar plot for n = 1 is given in 
Fig. 1^. When solving Eq. ( ^.58|) it turns out that there can be an imaginary part of 



on-shell for certain values of Pz- Let us take the R branch as an example. For Pz^ M. the 
hole branch can exist if Esn-i + Pz^n-i does not grow with p^, and that is possible since 



19 



s{E,p1 + eB{2n — 1)) can be negative when the logarithm in Eq. ( |5.54|) dominates, i.e. 



E > ypl + eB{2n — 1). The larger pz is, the closer E is to the above square root. This is 
why the hole branch goes exponentially close to the light cone for large Pz in the standard 
analysis fl^. Here we also have a factor {Esn — Pzi^n) in the right dispersion relation in 



Eq. ( |5.58|) . In Sn and the argument in the logarithms contains E — \jp^ + eB{2n + 1). It 
is then clear that the energy can go below the n-th light cone when pz is large enough and 
the logarithm picks up an imaginary part. The whole approximation breaks down at these 
points since the excitations are no longer stable. The value of the spectral function does 
not satisfy Eq. ( p.61| ) when this happens and it is not clear to what extent it is meaningful 
to continue beyond these points, or even to go close to them. 

6 Energy in the lowest Landau level — Anomalous 
magnetic moment 

The anomalous magnetic moment of electrons is defined from the energy spectrum of a 
particle nearly at rest in a weak field. Usually, it is calculated from the triangle diagram 
(see Fig. using plane waves as external states, while it really should be done in the 
proper Landau levels as explained in Section ^. Furthermore, not only is the shift of one 
level required but also the energy gap between different levels generated by the field. In 
vacuum the anomalous energy shift is independent of the Landau level for low enough fields 
(see e.g. We have not yet calculated the thermal shift for an arbitrary Landau level, 

but there are some interesting features already in the lowest Landau level. In vacuum one 
can calculate the anomalous magnetic moment either from a triangle diagram, using plane 
waves as external states, or extract it from the linear term in the exact expression with the 



electrons in the general Landau levels [^. The result is the same. At finite temperature 
the issue is more complicated because of the IR sensitivity and it really makes a difference 
if the external states are in a proper Landau level, or plane-wave states. A related problem 



was discovered in |^5|, ^ , where a real-time method was used for the triangle diagram, but 
we later found that the limit procedure of taking the external photon momentum to zero 
is not unique and thus the result is not well defined. We shall exemplify this phenomenon 
here by calculating the anomalous magnetic moment both from the expressions in Section 
^ and from the triangle diagram (Section in the static limit. 
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7 The Self-energy method 



Comparing with the energy shift for a charge —e, particle with spin (s) magnetic moment 
/X = —e/2mgs in a magnetic field B = —Be^, AE = — /x • B, we can write for an electron 
in the lowest Landau level ( here (s^) = 1/2): 

AE = (vl>i;o,o,o|S|vl>i;o,o,o) = Am^^^ - + 0{eBf , (7.62) 

where Am^'^ is the thermal mass, and 5g/2 is the anomaly of the magnetic moment. To 
linear order in the magnetic field, Eq. ( p.l3|) gives the very well-known result |12 



The total energy remains positive even for very large B when the non-linear terms dominate 
[ P^j . Let us now consider the thermal photon contribution using Schwinger's proper-time 
method. When the operators in Eq. ( |4.47D act on a Landau level, we can use the properties 
(see Appendix 0) 

^3*Sfo!o,o = ^Sfo!o,o, (7.64) 
]^±*!fo!o,o = 0, 

to compute the expectation value of T,^ in the lowest Landau level. In the limit of vanishing 
Pz we obtain to linear order in eB 

AB>. = 0) = ,n(-^ + ^— ^] . (7.65) 

which agrees with standard results [29-32], and also with the Furry picture in Eq. ( |4.42| ) 
and the limit ^ in Eq. ( [4.43| ). We also notice that the main contribution comes from 
the hard part of the loop integral (/c ^ T) and there is no particular IR sensitivity. 
For thermal electrons we similarly obtain in the limit of vanishing B and pz'- 

A£;f*_(B = 0,p, = 0) = Aml'X- 

= ^jyk,,e(kl-r.^)Mk,,)4W^^^^^ , (7.66) 
which agrees with the corresponding limit of Eq. ( [4.39| ). The thermal-electron contribution 



to linear order in the magnetic field is more involved. First, there is a 5 dependence in the 
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0-function in Eq. ( |4.51| ), which after expansion generates a derivative of the distribution 
function. Furthermore, one cannot expand the integrand in B naively, because that would 
generate IR-divergent integrals. Another delicate integration by parts of the logarithm in 
Eq. ( |4.51| ) is needed to give a finite integral after the expansion. There is thus an ambiguity 
in the way of writing the self-energy, depending on which terms are integrated by parts 
to produce derivatives of the distribution function. We have chosen to perform further 
integrations by parts to remove the derivative on the distribution function whenever this 
does not produce any IR divergences. After these manipulations we find, up to linear order 
in B: 

A J^p,^^ A (3,f, eB a duj 



2m Svr Jm — vn? 

' 2^"^ + 2muj — w? m \ „, . , ^ dftioj) 

Jf[uj) + 2m- 



rn? {ijj + m) J duj 



+ 



eB a f°° duj 2u'^ - Sm^u - 2rn? , , . 
2m 37r Jm ^/uj^ — m^ m'^[uj + m) 

In order to obtain this result using the Furry picture propagator, we must take the limit 
— ^ in Eq. (|4.41| ). Obviously the expansion will be in powers of Pz/iko — m). This 



produces IR singularities, if we are not careful. For the antiparticle part —ko = uj > 0, 
so here it is straightforward to perform the expansion. For the particle part we know, 
from above and by comparison with results using the imaginary-time formalism, that it 
is likely that the finite result should contain derivatives of the distribution function. Let 
us therefore introduce the factor of convergence {u — mY, where u is assumed to be so 
large that we may perform the expansion to get the contribution for vanishing p^. We may 
then isolate the term that will become divergent as z/ — > 0. Integrating this term by parts, 
the out-integrated term vanishes for z/ > 1. In the end we may consider the analytical 
continuation i/ — 0, and arrive at Eq. (|7.67|) . 
In the high-temperature limit we obtain 

Again, it agrees with standard results for the high-temperature limit and it is dominated by 
a hard thermal loop. The difference with Eq. ( |7.65| ) is that Eq. ( [7.68| ) is only approximate 
in the high-temperature limit and that sub-leading terms are IR-sensitive. This is why we 
had to be so careful with the expansion in B. 
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Figure 3: The triangle diagram. 



7.1 The triangle diagram 

The traditional way of computing the anomalous magnetic moment is from the triangle 
diagram in Fig. |^. In vacuum, the full vertex, sandwiched between plane-wave states, can 
be written as 



Up'T f,{p' , p)up = Upf 



2m 



Up 



Up' 



Up , (7.69) 



where the Gordon decomposition was used in the last equality. The anomalous magnetic 
moment is defined by Sg/2 = -^2(0) and it can be extracted from the term linear in in 
the limit q —>■ 0. When using the Gordon decomposition, it should be kept in mind that 
the external states are on the B = mass shell. The use of such states is questionable 
from the point of view of perturbation theory, as discussed in the beginning of Section ^. 

At finite temperature there are two formalisms for calculating the triangle diagram, 
the imaginary- and real-time formalisms (ITF and RTF). In the ITF we can put q^, = 
from the outset but we are left with an analytic continuation in the po variable. On the 
other hand, in the RTF (we shall use Thermo Field Dynamics (TFD) as a RTF) we get 
the result for real po, but it is more tricky to put q^ = before doing the loop integration, 
because there are potentially ill-defined products of distributions. There has recently been 
much interest in two issues that are of importance here. First, the relation between ITF 
and TFD has been clarified. It was first found that the naive diagrams in TFD (the ones 
with only physical fields on the external legs) did not give the same result as the ITF 
in general p9| , The difference was shown to be related to the different kinds of 

analytic continuations that are possible in the ITF, and different time-ordered, retarded 
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and advanced Green functions Q It is now clear that a certain combination of 

TFD diagrams gives the same resuh as the usual analytic continuation {uj —>■ pq + ie) in 

ITF, which is the retarded Green function. Secondly, the Lorentz invariance is broken 

at finite temperature so the vertex function depends on go and q independently. We can 

thus take the Qfj. ^ limit with either > or < (time-like and space-like limits) 

or generally with any fixed ratio |q| /go. It has been known for some time that e.g. the 

one-loop self-energy in the 0^-theory does not have a unique limit when g^ — ^ p5| , |46[ . 

In principle it would be interesting to discuss both limits, but the external mass shell 

conditions on p and p' put the constraint on g that 

1 
2 



Poqo - pq 



q 







(7.70) 



which prevents the use of go and q as independent variables. We shall therefore concentrate 
on the case with g = for both the ITF and RTF. 

In the ITF we can simply put g to zero from the outset. We use standard ITF rules in 
Euclidean space and write the expression for the vertex correction as 

2 



—I 



{k — pY + it \k'^ — m? + ie 
The simplified tensor structure in the numerator is 



T,{k) 



(7.71) 



(7.72) 



The last term can be replaced by 4(2m — k^jki, since we are only interested in terms 
proportional to pi (and not to 7i) and we neglect higher-order terms in pi in the non- 
relativistic limit. Then also 70 can be replaced by U. The double pole in k"^ 
conveniently be represented by a derivative. 



m can 



d 



dM^ k'^-M^ + ie 



m? + ieY 



(7.73) 



KP=m? 

After an analytic continuation we find as usual two contributions, one from thermal photons 
and one from thermal electrons. The photon contribution agrees with Eq. ( [7.65| ) so we leave 
that aside. From the thermal electrons we obtain 

5g%^^. a r-o _ ^ r 2 



-IT / ^^(^ 

OTT Jm 



m 



2\3/2 



2m — UJ 



duj muj{m — uj)'' 



+ 



dfp{uj) 2m + UJ 
duj muj{m + uj)' 



(7.74) 
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We would like to remark here that special care has to be taken when splitting into matter 
and vacuum contributions, when the Feynman diagram considered contains fermion as well 
as boson propagators. In the so-called finite-density contribution (see e.g. Eq. (3.71) in 



1^), there is a term 0(po ~ ^) that could have been mistaken for 0(/i — tu). This term 
has to be added to the vacuum contribution in order to give the correct Feynman pole 
prescription for the vacuum part of the fermion propagator. The leading high-temperature 



contribution comes from the highest power of u in the integral and it is given by 



9 \m 



5 



(7.75) 



which agrees with the result in Eq. (|7.68|) . The sub- leading terms, which are sensitive to 
the IR part of the loop integral, are on the other hand completely different. 

In TFD there are several different vertex functions depending on how the external legs 
are chosen to be particle or thermal ghost lines and we call them Aabc {a-ih^c = 1,2). 
According to [^, ^ Q it is necessary to sum over 1 and 2 with certain weights for all 



external points except one. The external point with only a 1-field is the one with the latest 
time argument. In our case we consider an incoming electron going through the heat-bath 
and interacting with an external magnetic field. The incoming point c in Fig. ^ must then 
have earlier time than point a. Also, the external field at point b does not infiuence the 
electron at later times than the time at point a. We are thus naturally led to consider the 



combination with a as the latest time. According to |^ the retarded Green function with 



the outgoing electron at the latest time and g = is given by 

A = (Am + Ai2i) + e't'P^^/\A^,2 + A122) • (7.76) 



In Eq. (|7.76|) there are several terms with overlapping 5-functions for = but they 



actually add up to give something finite. In particular we see that the incoming photon 
line at point b should be summed over 1 and 2 with the consequence that there is no 
product of 5-functions from the two fermion propagators. The 2x2 structure of the 
matrix propagator from c to a becomes 

/ G\k) - smH[G\k) - G*\k)] sin ^ cos ^[{G\k) -G*\k)] 

y - sin ^ cos ^[{G\k) -G*\k)] G*"^ (k) - sin^ ^[G^{k) - G*\k)] 

where G{k) = {k"^ — m'^ + ie)~^ andsin^i? = friko). The combination G'^ (k) — G*'^ (k) turns 



T,ik) 



(7.77) 
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out to be related to a derivative of a (5-function, 

6{k^ - M^) , (7.78) 



G^{k)-G*\k) = 2TTi ^ 



which eventually leads to exactly the same expression as in the ITF. The remaining multiple 
5-functions do not contribute for on-shell external particles when = 0. The only thermal 
contributions are linear in the electron and photon distribution functions. The same form 
of regularization of the product between a propagator and a 5-function with the same 
argument was advocated in but here it followed simply from using the rules of f^ . 
It has been shown for the n-point function in a scalar 0^-theory that the ITF (and the 
mass-derivative formula) is recovered in the space-like limit but not in the time-like ||45|| . 
This suggests that putting = from start corresponds to the space-like limit and the 
time-like limit could very well give a different result. For a static external field, with space- 
dependent gauge fields, it seems most appropriate to take the space-like limit. Nevertheless, 
we do not find agreement between the triangle diagram Eq. ( [7.74| ) and 5g^+^- inferred from 
Eq. (|7.67|) , apart from the leading high-temperature terms. We surmise that the reason 



for this is that perturbation theory in B is not adequate, as discussed at the beginning of 
this section, and that the full Landau levels should be used as external states. 



8 The strong-field limit 



In the limit of extremely strong magnetic fields, the energies in all but the lowest Landau 
level are approximately proportional to V eB. This implies that intermediate states of 
higher Landau levels are suppressed. It is thus very convenient to use the Furry picture 
propagator with its explicit spectral decomposition in this case. One may then neglect 
the contributions from all but the lowest Landau level, i.e. the sum over n is reduced to 



n = 0. In Section |8.1| we calculate the self-energy in this approximation. In Section ^.2| we 
explicitly verify that in this limit of strong magnetic fields the self-energy is independent 
of the gauge-fixing parameter in the photon propagator. 
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8.1 The self-energy in strong fields 



In it was shown that the approximation of only considering intermediate electrons in 



the lowest Landau level gives the same result as the high-field limit of Eq. ( p.l3| ), i.e. 

AB.„-|^f (lnH^|)%0(l„^) , eB»,n\A . (8.79) 
i^o 47r V / V / 

When the contribution from thermal electrons is considered, it is perfectly clear from the 
suppression by the distribution functions that the lowest Landau level is dominating for 
eB ^ m^,p^, /i^, T^. After reducing the sum to n = in Eq. (|2.14| ) and performing the 
integrals as before, we may now also perform the remaining Gaussian integrals over x,x' . 
The result reads 



Let us now use the 5-function to integrate over kz- Then we perform the angular integral 
in k_L and substitute u = k'j_. We find 



Bo 2ir i-oo ,/H _ ,„2 



POD 





exp(— -u) 



u - [{ko - EoY - {Vz ± V^o - m^? + ^e]/2eB 

We may here identify exponential integrals, and use their asymptotic expansions as given 
in Appendix Q The leading strong-field behaviour is thus 



vr J-oo ^/fcg - m? \m\ko-Eo 



m r°° duj „ , ^ , ^ 

Notice here the field-independent (in this limit) imaginary part, which is exponentially 
suppressed for weaker fields. We have also investigated the self-energy for a positron in 
the lowest Landau level. The result is the same as in the electron case above, but fp{uj) 
is replaced by fp{oj) in the imaginary part on the right-hand side of Eq. ( ^.82|) . The 
physical process which is responsible for the occurrence of the imaginary part of AE'^f'^. 
is the reaction e~ ^— 7, which is possible in the presence of an external field that may 
absorb momentum to make the process energetically possible |5^. We have calculated in 
Appendix ^ the tree-level decay rate F, for an electron in a background of positrons. The 

27 



result is just as expected, T = 2ImA£'g|g_. The way we obtain the imaginary part in 
AE^+'^- in the strong-field limit with the +e prescription in Eq. ( p.81| ) is similar to the 
pole prescription of Landau in the theory of longitudinal plasma oscillations (see e.g. [pi]]). 

For the thermal-photon contribution we use instead the 5-function to integrate over k±, 
and find the result corresponding to Eq. ( p.81| ), 

AE^^ ^ p r ikMk.) : . (8.83) 

where k = |/co|. We may, to leading order, approximate the exponential function with 
unity. The poles of the integrand are outside the interval of integration in A;^, so the 
thermal photon contribution is real, and we may let e vanish. However, k^ + 2£'o^o + p1 
changes sign at k^ = —Eq ± m, so we must split the k^ integral in two parts. After 
performing the fcz-integral we thus obtain 



EonJEo-m j2Er.k- k^ -pI 

+ arctan 




k-Pz 



2Eok -k'^-pl 

dkofniko) 



X In 



Eq 271 [J^oo J-Eo+mj ,JW+2EQki+Pl 

e + E^ko + kJk^ + 2Eoko+pf 



k^ + Eoko - kJk^ + 2Eoko + pI 



(8.84) 

Notice that this is finite since there are cancellations between all of the ostensible IR 



divergences. For > 0, fco = is contained in the last integral in Eq. ( 8.84 ). To leading 
order the integrand is odd in ko for small ko, and a cancellation thus occurs. The next 
leading term produces a finite result when integrated over ko. For Pz = there is a 
cancellation between the leading terms from the first and last integrals in Eq. ( 8.84|) , and 
the result is still finite. 

We have plotted the different contributions to the self-energy in the high-field limit 
{eB ^ m^, pI, /i^, T^}: in Fig. ^ as a function of the magnetic field, in Fig. ^ as a function 
of the temperature, and in Fig. ^ as a function of the momentum in the 2;-direction parallel 
to the magnetic field. The self-energy is even in pz. In each case the self-energy is small 
compared to the electron rest mass and thus will only provide a small energy shift. In the 
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Figure 4: The electron self-energy AE = AE^vac + ^E^e+e- + high-Reld limit 

{eB 3> m^, pI, /i^, T^}, as a function of the magnetic field, for Pz/m = 1, iJ,/m = 1, and 
T/m = 1. 



high-field limit it is the vacuum contribution that dominates. Using Eq. 



we see that 



for the self-energy to be of the same order of magnitude as the electron mass we must have 
eB/m? ^ 10^''. This corresponds to the mind-bogglingly large field oi B c^ii 10^^ T. 



8.2 The gauge-fixing dependence 

We formally show in Appendix that the self-energy is independent of the gauge-fixing 
parameter ^ on the tree-level mass shell. In the limit of extremely strong fields we may 
explicitly verify this conclusion. We shall here consider the different contributions to the 
self-energy appearing for non- vanishing ^ in Eq. ( |2.7| ). To regularize the vacuum contribu- 
tion we assume the external electron to be off-shell 



7-12 2 2 

h — m — 



A 



^5) 



and we shall find that the gauge- dependent part is proportional to A, as shown formally 
in 



52| . Factorizing out ^, the vacuum contribution reads 



2eB 



dkn dky 



due 



X 



i2Tif2m{E + pz}J-oo JO 
[{ko + kz)A^ + {kl - - kl + ie)[-2E^ - {E^ - m^)kz + 
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Figure 5: The electron self-energy AE = AEyg^ + AE^f^_ + AE^ in the high-held limit 
{eB S> m^, p1, /i^, r^}, as a function of the temperature, for eB/m^ = 8, Pz/m = 1, and 
ji/m, = 1. 
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Figure 6: The electron self-energy AE = A£^vac + ^E'^+e- + ^^l^ ^'^ the high-held limit 
{eB S> m^, pI, ij?, T^} as a function of the momentum parallel to the magnetic Geld, for 
eB/m^ = 8, ^/m, = 1, and T/m = 1. Each contribution is even in p^. 
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+ {E^ + m2)/co + 2p,{E{ko - h) - + pi} + pl{2E + h - h 

1 

^ (ki -m?-kl + ie)[{ko - Ef - {k, - p,y - 2eBu + ie] ' 
Let us now use the Feynman parametrization 



in the first term on the right-hand side of Eq. ( |8.86| ). We may then integrate over ko and 
kz, to obtain 



A a r°° r 1 

AE,ac,e(5 ^ oo) ~ -— / rfwe-" -+ 

m Air Jo lu — le 



2eBA / ds 



1 



lo [m2(l-s)2- As(l-s) + 2e5Ms-ie]2j ' ' 

Notice that we cannot naively let A = here, since the result would be an ill-defined 
product of A and a logarithmic divergence. The u integral is dominated by small u. The s 
integral would be dominated by s ~ 1. Inserting s = 1, except in the terms (s — 1), we may 
easily perform the s integral. The result will contain simple polynomials and logarithms. 
Expanding the logarithms for ura^/eB ^ 1 to leading order, several cancellations will 
occur. The result reads 

AE^acf(5 ^ oo) ~ — / du — — . (8.89) 

^'^^ ^ 477 m(2m2 - A) io u+ (m^ - A) /2eB - ie ^ ' 

We may here identify an exponential integral, and use its asymptotic expansion to get 

AB,^,(B-oo).---p-^ln(-^) . (8.90) 

This is obviously vanishing on the tree-level mass shell A = 0. However, notice that if we 
would like to solve the dispersion relation self-consistent ly in this limit, the gauge-fixing de- 
pendence is only logarithmically suppressed compared to the result in the Feynman gauge, 
proportional to \\.\i{2eB / rn^)]^ according to Eq. ( ^.79|) . When considering the contribution 
of thermal electrons to the ^ dependence, the integrand in correspondence to Eq. ( p.86| ) is 
proportional to 

S^kl -m"- kl){{ko + kz)A? + {kl -m"- kl)[-2El + ...]} . (8.91) 

This is vanishing on the tree-level mass shell A = 0. Also for thermal photons we may 
immediately use the on-shell energy. The contribution to the self-energy multiplied by ^ 
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then becomes 

AE^.c^-eB^ due-"" dkodKko^'ikl - kl ~ 2eBu)fB{h) , (8.92) 
''^ 7r m JO J~oo 

where the prime on the 5 function denotes a derivative with respect to its argument. Since 
the photon is its own antiparticle, /b(/co) has to be even in k^. Then AE^ ^ is vanishing 
due to the antisymmetric integration in k^. 

9 Discussion and final remarks 

In this paper we have presented a detailed account on the fermion self-energy when both 
an external magnetic field and the presence of a heat bath have to be taken into account. 
In the context of astrophysics, many QED processes in the presence of strong magnetic 
fields have been studied, however mostly without the presence of a heat bath. 

As was mentioned in Appendix O it has been argued in the literature that one must 



use a Dirac spinor basis which diagonalizes in the Sokolov-Ternov spin operator [53 



in order to consistently treat unstable excited relativistic Landau levels in perturbation 



theory [^, We have verified that such a basis can be constructed from the Dirac 
spinor basis defined in Appendix ^ We have also seen that when a heat bath is present 
the situation is more complex and that in general such a basis does not diagonalize the 
self-energy except in the case of zero momentum parallel to the external magnetic field. 

The effect of thermal quasi-particles on the emission of neutrinos of a very hot or dense 
star, but without the presence of a magnetic field, has been studied in [0. It was found 
that under certain physical conditions the electron-plasmino annihilation may exceed the 
electron-positron process into neutrinos. In principle it is possible to carry out a similar 
analysis in the presence of a magnetic field. As we have seen in Section ^ the presence 
of a magnetic field leads to a pz asymmetric quasi-particle spectrum. We have suggested 
that this asymmetry may play a role also in dynamics, leading to the observed high space 



velocities of pulsars |^, |34| . 

We have also carried out a detailed analysis of the anomaly of the magnetic moment 
5g/2 when a thermal heat bath of photons, electrons and positrons is present. Previous 



considerations ^ of this problem have led to ill-defined IR behaviour. With ther- 
mal photons, electrons and positrons present, we obtain the following contribution to the 
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anomaly: 



2an T2 
9 



duj 



X 



2uj^ + 2mco' — m? 



m 



a 



duo 



{uj + m) 
2uj'^ — Sm^u — 2m^ 



2m 



(9.93) 



■m? ^^{uj + m) 

To the best of our knowledge this is the first IR well-defined expression for the anomalous 
magnetic moment when a thermal environment of photons and electrons is present. At 
zero temperature the tu-integral in Eq. ( |9.93| ) can be carried out explicitly, and we find 



^9, 



e+e~ 



a 
3^ 



2 + 



m 



m 



m + /i 



mj 



2m 



m^ 



(9.94) 



for > m. In the limit of large fi, this agrees with the result in pH], as could be expected 



since the IR sensitivity is sub-dominant. Thus the conclusion in |]35[ that g can become 
considerably smaller than 2 for high densities remains valid. Using a bold extrapolation to 
large corrections we find that g + 6g c:^ for fi 35 m^, which should be compared with 
a typical chemical potential fi ^ 300 me inside a neutron star. Even if the approximations 
are not valid for such a large /i there could still be important corrections to the synchrotron 
radiation from the surface of a neutron star, and thus to the estimation of the S-field from 
observations. 

We notice also in Eq. ( |9.94| ) that Sg seems to become very large for small densities, e.g. 
when fi —>■ m. This is an artefact of the expansion in B. As discussed in Section ^ we do not 
expect an expansion in B to be universally possible, but some sort of de Haas-van Alphen 
oscillations should show up in certain limits. For T = and B fixed we can take the limit 
H ^ m directly in Eq. ( |2.14| ) and it corresponds actually to a strong-field calculation since 
only the lowest Landau level is inside the Fermi sea. The leading contribution can then be 
extracted from Eq. ( ^.81| ) using Eq. ( |E.31| ) and we obtain 



AE 



a 



TT 



m^ 



In 



2eB 



^2 _ jyfi 



C 



o 



m 



eB 



(9.95) 



valid for m < < y/m? + eB. Our conclusion is that the energy shift is perfectly finite 
when /i — i> m (in fact, it goes to zero), but since AE does not admit a power series 
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expansion in B the standard definition of a magnetic moment (see Eq. ( [7.62| )) does not 
make sense. 

In a strong magnetic field we have, finally, observed the presence of an imaginary part 
in the self-energy which becomes field-independent if the external magnetic field is large 
enough. Formally, this imaginary part appears in a way similar to the pole description of 
Landau in the theory of longitudinal plasma oscillations ||5l|]. This imaginary contribution 
is therefore a relativistic counterpart of Landau damping and is physically due to the 
possibility of e~^e~ annihilation into a single photon in a magnetic field ISC 
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APPENDICES 
A External-field propagator 



In this appendix we summarize, for the convenience of the reader, the relevant expressions 
used in the main text for a constant magnetic field B in the negative 2;— direction in the 
gauge = (0, 0, Bx, 0). We use the 7-matrices in the chiral representation, 

^"=(-11 j ' ^ =[-a^ j • ^^-1) 
We seek solutions to the Dirac equation for a fermion in an external field 

(^^-m)^W(x,t) = , (A.2) 
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where ilf)= 'y^{id^ + eA^) for a particle of charge — e. With our choice of gauge it follows 
that 



vl/W(x,t) 



1 



e^^[±i{-Et+pyy+p,z)]^f\x) , 



(A.3) 



2'kV2E 

where k denotes Py, Pz, and any other quantum-number necessary to specify the wave 
functions. For ease of notation we shall write iD^ — 11^ = {E, —pz, — H^) when acting on 
the above wave functions. Acting with (fl+m) on Eq. ( |A.2| ), using [n^;, Uy] = ieF^y = ieB, 
we find 

(E^ -m^ -pI-UI + eBa^y)^ = 



(A.4) 



where a^y = i[7a;, 7j/]/2 = diag [cr^, cr^]. Let us now introduce the functions 



TT / 



exp 



2 V eBj 



where Hn is the Hermite polynomial given by the Rodrigues formula as 

HJx) = 



1„2 C?" —Ix^ 



-l)"e2" — e 2 



and we define I^i-p^^x) = 0. The functions In-p^^x) are normalized according to 



(A.5) 



(A.6) 



(A.7) 



for ra, n' = 0, 1, 2, . . . . Defining S,± = T^xT- i^y, we have II^^ = + eB, and they act on 
the /„ functions according to 



i+In-i-pAx) = iV2eBnIr, 



IPy 



X 



^-In-py{x) = -W2eBnIn^i.py{x) 



(A.8) 



It readily follows that 



- eBaz) diag[4;p^(x),/„_i;p^(x)] = 2eBndiag[In;py{x), In-i;py{x)] 



We may therefore write $ in Eq. ( [A.3|) in the form 



*cSp„p.(^) = diag[/n;p,(x), /„;p,(x), J„_i;p^(x)]M^^j 



(±) 



:Py,Pz ' 



(A.9) 



(A.IO) 



where m^"^^ is a Dirac spinor independent of x^, and ( = ±1 denotes a polarization index. 
The energy eigenvalues are given by the relativistic Landau levels 



E = EniPz) = ^rn? + pl + 2eBn 



(A.ll) 
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In the lowest Landau level (n = 0) 
Uk possible, corresponding to C = 1 



implies that there is only one choice of 
For the higher Landau levels there is a twofold 



degeneracy, since two linearly independent u^^^^p^^^^ may be found, and there is thus an 
ambiguity in the choice of wave functions. 



In there is a discussion about a choice of the u^^] „ „ which diagonalizes the self- 



energy corrections in vacuum. We have not found any simple generalization of that basis at 
finite temperature, as stated in Appendix 0. On the other hand, in this paper we need no 
such basis since we always deal either with the expectation value of the self-energy in the 
non-degenerate lowest Landau level or with the full 4x4 matrix. For any higher Landau 
levels it is, of course, mandatory to diagonalize the matrix including corrections since the 



perturbation theory is degenerate, as emphasized e.g. in [14] 



We shall here use the choice of wave functions given in [|T2 



{En+Pz)L,r,AX] 
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(A.12) 
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y/En - Pz 
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' x\ 



(A.13) 
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(A.14) 



/ iy/2eBn In 
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'-En+Pz)In-l;-py[X 



-mln-i-py{x) j 



(A.15) 



It can be shown that the collection of all \E''s forms a complete orthonormal set. The wave 
functions are normalized according to 



/ d=^a:^(;!;^,,^(a:)v^(^3,^, (x) = <5c,C' Sx,x' S{py - p'^) S{pz 



Pz 



(A.16) 



where A = ±. The propagator in the external field iS{x',x) = «5'vac{x',x) + iS^'^{x',x), 
as given in Eqs. (p.3|), ( |2.4| ), is then explicitly written (see e.g. [1^ for the vacuum part. 
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and for the additional term when we have fermions according to some one-particle 
distribution): 



S{x', x] 



ab 



E / ^^^^ exp[-tko{t' -t) + tky{y' -y) + z^z' - z)] 



n=0 
X 



(2vr)= 



1 



kQ — k1—iv? — 2eBn + ie 
^Sabi^'i kQ, ky, kz'i X , x) , 



+ 2m6{k^-ki-m^-2eBn)fFiko 



where 



fpiko) = e{ko)fUko) + Q{-h)fFi-ko 



(A.17) 



(A.18) 



The matrix S{n; ko, ky, k^, x', x) entering above is in the chiral representation of 7-matrices 
explicitly written 



S ("Ti; ko, ky, kz 








- {ko + k z)In,r. 

mln-i ,n-i iV2eBnIn-i,, 
-{ko-kz)In,n i\/2eBnIn ,n— 1 
\ -W2eBnIn-i,n -{ko + kz)In-i,n-i 



where we have used the short-hand notation 



B The Sokolov— Ternov spin operator 



-iV2eBnIn,n-i \ 

■{ko — kz)In-l,n-l 



mln-l,n-l ) 



(A.19) 



(A.20) 



Let Hd = a • n + /3m, where a* = 7*^7*, (3 = 7*^, be the Landau-Dirac Hamiltonian in an 
external magnetic field. Let the magnetic field be parallel to the z-axis. The z-component 
of the Sokolov-Ternov ||^ spin operator is fiz = o'zijn + (For a recent discussion of 
this operator in the context of calculations of the electron anomalous magnetic moment, 
|5^.) It follows that fiz is conserved, i.e. it commutes with Hd. It is easy to construct 



see 



linear combinations of the solutions in Eq. ( |A.3| ), which also diagonalize fiz- Suppressing 
all the quantum numbers of the solutions in Eq. ( |A.3|) except the polarization = ±1, 
solutions to the Dirac equation ( [A.2|) that diagonalize fiz can be written in the form 



(±) 



{B.21] 
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where, if n 7^ 0, 



where p„ = \/2eBn. The states in Eq. (|B.21| ) have the property A^^/i^'' = ±Jrn^~+^'^^^^ 



The lowest Landau level is diagonal in /t^ with eigenvalue m. It can be verified that the 
propagator S{x',x), when calculated in the basis in Eq. ( p.21|) , still is the same as the 
one given in Appendix 0. The basis in Eq. (|B.21|) furthermore has the property that the 
self-energy is diagonal in the vacuum sector to 0{a). By using the fact that the 

on-shell Sokolov-Ternov operator can be written in the form fiz = (^zil^P^ — IzPz), this 
follows easily from the fact that [fiz, P^o] = 0, where Sq is the vacuum self-energy operator 
[^]. In the presence of a thermal heat bath the self-energy operator defined as the 

sum of Eq. ( p.32| ) and Eq. ( p.34| ), commutes with fiz only if = 0. It has been observed 
|5^ that in the vacuum sector one must use the basis where fiz is diagonal in order to 
introduce consistently in perturbation theory, a resonance width into the QED cyclotron 
scattering amplitudes. In a thermal environment the situation is more complex and further 
studies are required. 



C The tadpole 

Also the tadpole could possibly contribute to the one-loop self-energy. In configuration 
space the tadpole is proportional to 

r{x) =tT[eYiS{x,x)] , (C.23) 

where the trace is over spinor indices. Introducing an e > in the argument of \1/ : t' — 
t + e, we may remove the time-ordering. Then we let e vanish and immediately find 

r{x) = -e(^{x)Y'^{x)) , (C.24) 

i.e. the expectation value of the electromagnetic current. We could instead choose to change 
the argument of \1/ : t — > t+e, in order to remove the time ordering, and then let e vanish. In 
this case we need to use the equal time anticommutation relation {\l/a(x', t), \I'|^(x, = 
SabS^i'x' — x) to reverse the order of \E' and Since 7^^ is traceless, we again arrive at 
Eq. (|C.24|) . Let us now separate the current in its vacuum and thermal contributions 
= + J^^^. Trivially J^^^ vanishes after renormalization, and so does J^^ , for 
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j = 1,2,3 in our static model. What is left is the charge density Jg_^. But in order 
for our static model to be valid there must be no such charge asymmetry on average. If 
there is a finite chemical potential for electrons there must be a compensating (static) 
background charge to make the average charge density vanish everywhere. Actually, it is 
easily shown with explicit calculations for the particular system we consider that the only 
possible non- vanishing part of J'^, as defined in Eq. (|C.23|) , is J°^. Using the electron 



propagator in Eq. ( |A.17| ) and its thermal generalization, it follows that and vanish 
when performing the trace. Also and J^^^ vanish due to the antisymmetric integration 
in ky and ko, respectively. In the case of vanishing chemical potential, friko) is even, and 
then also Jg ^ vanishes due to the antisymmetric integration. In the case of general fi we 
find 

4,, = 7^ E / dk4f^[E4k.)] - /^[E„(A;,)]}(2 - <5„,o) , (C.25) 

i.e. the charge density. Notice the factor (2 — (5„ q); which originates in a term In,n+In~i,n-i, 
and shows the twofold degeneracy in all but the lowest Landau levels. 

D Gauge independence 

We have used the Feynman gauge throughout the article and now we shall show that the 
result is gauge-independent on-shell, i.e. (^'kISI^'k') is independent of the gauge-fixing 
parameter ^ when and \E'k/ are solutions to the Dirac equation. The only possible ^ 
dependence in a general covariant gauge comes from the photon propagator. In vacuum, 
it would give a term of the form 

e J d^x' ct^x^ ^{x')ie^-i^'[d^dlb{x' - x)]S{x, x')Y^k'{x) , (D.26) 

where 

After one partial integration in x or x' the Dirac equation can be used together with 
{ip— m)S{x,x') = 6{x — x') to show the Eq. ( p.26| ) is zero. At finite temperature the 



gauge-dependent part of the photon propagator is again given by a total derivative 
and the same proof goes through. 
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E Exponential integrals 



In [515 1 we have the following definitions of the exponential integrals 

TOO 

Ei{z) = / dt — ,|argz|<7r , (E.28) 



oo 



Ei(x) = -V j dt^ , xeJR . (E.29) 
Furthermore they satisfy the interrelation 

Ei(-x±i£) = -Ei(x) T27r ,x>0 , (E.30) 

and have the series expansions 

oo / -\\n^n 

Ei{z) = -C-ln;z-V^ ,|arg2|<7r , (E.31) 

n=i ^'T^- 

oo n 

Eifx) = C + lnx + y — r ,a;elR , (E.32) 
where C = 0.5772156649 ... is Euler's constant. 

F Electron-positron annihilation 

We shall here consider the process e^e+ ^ 7, where the positron comes from the heat and 
charge bath. A large magnetic field may absorb momentum to make the process energeti- 
cally possible. Consider an electron in the state described by a positron in the state 
described by '^[r\ and a photon with momentum and polarization e^{X, q). To the low- 
est order in perturbation theory (i.e. perturbatively in the quantum electromagnetic field, 
the static uniform magnetic field is treated exactly), we find the corresponding transition 
matrix element 

^Tp,«;,,A = -le j d''x^^-\x) ^^¥;-\x)el{q, \)e''^-^ . (F.33) 
Let us now form |Tp K;g,Ap- Sum over all polarizations of the outgoing photon 

and integrate over the photon momentum 

/ ^ = / [dH 2n6iql ' q') ©(^o) • (F.35) 



go=|q 
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Sum over all incoming positrons with distribution fp{Ef^). Due to the normalization of 
our wave functions in Eq. ( |A.16| ) we must divide by / dydz, for the norm of the decaying 
electron state. Dividing by the infinite time elapsed between the final and initial state / dt, 
we find the decay rate 



/ dydz dt 

X /[A]2vr5(g2)e(go)e-^-(^-'-^)vl/W,^,,^(x') . (F.36) 

By comparing with Eq. (|2.4]) we see that /f (-^k)^! \^) is exactly the thermal 

positron contribution to the Dirac fermion propagator. Due to energy conservation (as 
follows when performing the integration over t) go < will not contribute, so we may drop 
0(go)- The decay rate is thus exactly equal to the imaginary part of the electron self- 
energy, with thermal positrons in the intermediate states, obtained through 1/ (g^ + ie) = 
V^i/q^) — i'n:5{q^) in the photon propagator. This splitting into principal and imaginary 
part is equivalent to what was done in the high- field limit in Section |8]^, using the definition 
and interrelation of exponential integrals as given in Appendix |^. 
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